Quantum mechanical photon-count formula derived by entangled state representation 
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By introducing the thermo entangled state representation, we derived four new photocount dis- 
tribution formulas for a given density operator of light field. It is shown that these new formulas, 
which is convenient to calculate the photocount , can be expressed as such integrations over Laguree- 
Gaussian function with characteristic function, Wigner function, Q-function, and P-function, respec- 
tively. 

PACS numbers: 



In quantum optics photon counting is important for 
judging the nonclassical features of hght field, most mea- 
surements of the electromagnetic field are based on the 
absorption of photons via the photoelectric effect. This 
is true not only for used insofar as photodiodes, pho- 
tomultipliers, etc., but also for such homely devices as 
the photographic plate and the eye. So the problem of 
photo-electric detection attracts an increasing attention 
of many physicists and scientists. Expressions for the de- 
tection probability have been presented in many works 
[l|, [2| • The quantum mechanical photon counting distri- 
bution formula was first derived by Kelley and Kleiner Q . 
As shown in Refs. for the single radiation mode, the 
probability distribution p (to, T) of registering m photo- 
electrons in the time interval T is given by 
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where the state vector |^) is conjugate to the state \ri) , 



" is a displacement operator, and is 



a fictitious mode accompanying the real photon creation 
operator a* , |0, O) = |0) |0) , and |0) is annihilated by d 
with the relations [a, a^] — 1 and [a, a^] = 0. It is easily 
seen that \ri — 0) and |^ = 0) have the properties. 
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where C oc T is called the quantum efficiency (a mea- 
sure) of the detector, and : : denotes normal ordering, 
p is a single-mode density operator of the light field con- 
cerned. The aim of this Letter is to derive some other 
quantum mechanical photon-count formula by introduc- 
ing the thermal entangled state representation and con- 
vert the calculations of Wigner function (WF) and the 
characteristic function of density operator to an overlap 
between "two pure" states in a two-mode enlarged Fock 
space, so that it is convenient to calculate the photocount 
when a light field's density operator is given. In addition, 
this new method seems concise and easy to be accepted 
by readers. 

Recall that the thermal entangled state representa- 
tion (TESR) is constructed in the doubled Fock space 

0] based on Umezawa-Takahash thermo field dynam- 
ics (TFD) i.e.. 
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where n = n, and h denotes the number in the fictitious 
Hilbert space. 

According to the TFD and Eq.(l4|), we can reform the 
probability distribution p (m, T) as 
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where in the last step, we have used the operator identity: 
exp (Aatfl) =: exp [(e^ — l) a^a] : . Note that the den- 
sity operators p{a^ ,a) are defined in the real space which 
are commutative with operators {d\d) in the tilde space 

with \p) = p\I) , as well as (n| tj — Sn,i {n = h,l = 

I). By using a"' \l) = y/n/{l - my.\l - m) ,a^"'- \l) = 
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+ my./U |; + m) , Eq.® becomes 
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In order to derive four new formulas for p {m,T), we 
first bridge the relation between the characteristic func- 
tion (CF) and the entangled state representation {ri\. 
Similarly to Eqs.®, after using the TFD theory, the CF 

of density operator p, xs (A, A*) = tr ^pe'^°^~'*'*°^ , can 

be calculated as 
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which is the CF formula in thermo entangled state repre- 
sentation, with which the characteristic function of den- 
sity operator is simplified as an overlap between two 
"pure states" in enlarged Fock space, rather than using 
ensemble average in the system-mode space. Thus we 
can then simplify the calculation of xs (A, A*) by virtue 
of some important properties of the entangled state rep- 
resentation (77 1 . 

Using the expression of {ri\ in Fock space, i.e.. 
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where Hm n {S,* , £,) is the two- variable Hcrmitc polynomi- 
als [nl, [l3 , one finds 
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which leads to 
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where in the last step, we have used the formula |13l |. 
^ -j^H„,+i^„+i (x, y) 
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and the relation between two-variable Hermite polyno- 
mials and Laguree polynomials. 
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Further inserting the completeness relation of (77I ,i.e., 
J id^?7 \r\) (ryl = 1 (it can be proved by us- 
ing the normally ordered form of vacuum projector 
|0, 0) (0, 0| = : exp (— a^a — a^a) : and the technique of 
integration within an ordered product (IWOP) of opera- 
tors [l4l - [l6j ). into Eq.(l7]), we can rewrite it as 



p (TO, T) = 1 1 Q lA^) XS (A, A*) , 

(14) 

which is just a new relation about the CF and the photon- 
count distribution. When the characteristic function 
Xs(A, A*) of density operator for Wigner-Weyl form is 
known, the photocount distribution can be calculated by 
using Ea.dltl). 

For instance, we first consider the single-mode coherent 
states 1/3), whose CF reads 
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substituting it into Eg. p^ yields 
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where we use the limiting expression 
liuix^Q x"^ L„i{— \a\'^ / x) = ^ |q;|^™ and the follow- 
ing integrational formula (see Appendix B), 
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Eq. dTB)) is the Poisson distribution coinciding with the 
result in Refs. 

As another example, we consider the single-mode 
squeezed vacuum state, exp [r (^a^^ — a^) /2] |0) , whose 
CF reads 
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substituting Eq. ([T8| into (fT4|) . we have (Appendix C) 
^"sechrtanh'^r „ / V 
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which Pm {x) is the Legendre polynomial and Eq. p^ is 
a new result. 
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Next, we derive other three new formula. Notice that 
the characteristic function xs{^t^*) is related to the 
Wigner function, Q-function and P-representation by the 
following Fourier transforms. 



X5(A,A*) = 
Xs(A,A*) = e 
Xs (A, A*) 



^Aa* —A* a 



W{a) (fa, 



(20) 



e 2 



Q{a)(fa, (21) 
^P{a)cPa, (22) 



respectively, thus substituting Eas. (pn)) - (P^ into (HH we 
can directly obtain 
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where (a) = 2tr (pA (a, a*)) ,Q (a) = i (a| p |a), and 
the integrational formula ((T7)) is used. Eg s. ([251) - ([25)1 are 
the new formula for evaluating photon count distribu- 
tion. Therefore, once one of these distributions of p is 
known, the photocount distribution can be calculated by 
using Eq. (P5)) - (l25t . which involve the Wigner function, 
Q-function, and P-representation of p, respectively. To 
confirm their correctness, we still consider the coherent 
light field |/3) (/3| , its Wigner function, Q-function and 
P-function are given by W {a) — ^e"^'^""' , P (a) = 

5^"^^ {P — a), and Q [a] — ■^e^''^^"!^, respectively, then 
according to ([23)) - (l25]) and using p7)) and the above limit- 
ing expression lYm^^i:) L,n{— \oi\ /x) — ^ |q;| ™, one 
can draw the same result as Eg. dTHl) . 

At last, we should mention that using Eqs. 0-® 
it is shown that the Wigner function of a mixed state 
p, Wp (a, a* ) = 2tr ( A (a, a* ) p) where A (a, a* ) is the 
single- mode Wigner operator [13, 13 : whose explicit nor- 
mally ordered form is [l^ 

A (a, a*) = i : e-2(at-a*)(a-a) . ^ 1^ ^ga) (-l)°■^^ 

TT TT 

(26) 

which can also be converted to a overlap between two 
"pure state" in the enlarged Fock space. 



which is the Wigner function formula in thermo entan- 
gled state representation, with which the Wigner function 
of density operator is simplified as an overlap between 
two "pure states" in enlarged Fock space. Employing its 
completeness, i.e., J ^ |^) (^| = 1, one can derive these 
above new formula. In addition, the expression in Ea. (j27p 
can also examine the evolution of Wigner function of den- 
sity operator interacting with the environments (20| . 

In summary, based on Umezawa-Takahash thermo field 
dynamics theory, after introducing the thermo entan- 
gled state representation, we converted the calculation of 
CF to an overlap between two "pure states" in enlarged 
Fock space. Then we bridge the relation between the 
characteristic function and the photo-count distribution. 
Once the CF of density operator for Wigner- Weyl form 
is known, the photocount distribution can be calculated 
conveniently. Using the Fourier transform relation be- 
tween the CF and the distribution functions, we further 
derive other three new formula so as to be convenient 
for calculating photo-count distribution by using these 
formulas. 
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Appendix A: Derivation of sum-formula in 
Eq. (fT2)) 

Using the integration of two-variable Hermite polyno- 
mials. 
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Then making scale transform and using Eq.(Al) again, 
Eq.(A2) can be put into the following form 
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Appendix B: Derivation of integration- formula 
in Eq.dnD 
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Using Eq. (IT3l) and the generating function of the two- we find 
variable Hermite polynomials, 
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where we have used the formula 
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Using Eqs.(Bl) and ([T^ again, one can get the 
integration- formula in Ea. (|17p . 



In order to obtain Eq. (|19p . we first derive a new inte- 
gral formula, 
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Using Egs. p^ and (Bl), Eq.(Cl) can be put into the 
Appendix C: Derivation of the result in Ea.(|19p form 
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where in the last step, we used the formula [21| 
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whose convergent condition is Re{(zt f ±g) < 0, Re[(C2 - 
4/5)/(C±/±g)] <0. 



(C3) Expanding the exponential item involved in Eq.(C2), 

we see 
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where 
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Recalling that newly found expression of Lagendre poly- 
nomial (it is equivalence to the well-known Legendre 
polynomial's expression [22|). 



[m/2] 



1-4) ^Prn{x), (C6) 



the compact form for / is written as 
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which is a new integration formula. 
Substituting Eq.lHH]) into dH]) we have 

where /' shown in Eq.(C7) characteristic of 
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which leads to 
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which is the photon-count distribution of squeezed vac- 
uum state. 
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